REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  07 04-0  U  3 


The  public  reporting  burden  tor  this  collection  ot  intormation  is  estimated  to  averege  1  hour  per  response,  including  the  time  tor  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  meintaining  the  data  needed,  and  completing  end  reviewing  the  collection  ot  intormetion.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of 
intormation,  including  suggestions  tor  reducing  the  burden,  to  the  Department  ot  Defense,  Executive  Services  and  Communicetions  Directorate  (0704-0188)  Respondents  should  be  aware 
that  notwithstanding  any  other  provision  ot  law,  no  person  shall  be  subject  to  eny  penelty  tor  tailing  to  comply  with  e  collection  of  information  if  it  does  not  display  a  currently  valid  OMB 
control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ORGANIZATION. 

1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

08-10-2009  Journal  Article 

3.  DATES  COVERED  (From  -  To) 

4.  TITLE  AND  SUBTITLE 

Assimilation  of  HF  radar-dcrived  radtals  and  total  currents  in  the  Monterey 

Bay  area 

5a,  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

0601 153N 

6.  AUTHOR(S) 

IgorG.  Shulman,  Jeffrey  D.  Paduan 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

73-8404-07-5 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory 

Oceanography  Dtvtsion 

Stennis  Space  Center,  MS  39529-5004 

8.  PERFORMING  ORGANIZATIO  i 

REPORT  NUMBER 

NRL/J A/7330-06-7  30 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Offtcc  of  Naval  Research 

800  N.  Quincy  St. 

Arlington,  VA  22217-5660 

10.  SPONSOR/MONITOR’S  ACF  DNYM(S) 

ONR 

11.  SPONSOR/MONITOR  S  REP  )RT 

NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  is  unlimited. 


13  SUPPLEMENTARY  NOTES 


14  ABSTRACT 


20090814036 


Impact  of  I  IF  radar  surface-current  assimilation  on  ocean  circulation  model  predictions  in  the  Monterey  Bay  area  is  studied  and  evaluated  during  t  e  time  frame  of 
the  Autonomous  Ocean  Sampling  Network  (AOSN-II)  experiment  (August  -September  2003).  In  the  first  instance,  a  previously  described  method  ft  *  assimilation  of 
surface  current  data  is  applied  to  33-h  low-pass-filtered  data  and  a  non-tidal  version  of  the  circulation  model.  It  is  demonstrated  that  assimilation  ol  surface  velocity 
data  significantly  improves  the  surface  and  subsurface  correlation  of  model  currents  with  moored  current  observations.  These  results  from  the  \OSN-ll  period 
illustrate  that  surface-current  assimilation  is  beneficial  even  in  cases  for  which  very  high-resolution  (3  km)  atmospheric  forcing  is  utilized.  The  assin  ilation  approach 
is  also  tested  with  hourly,  unfiltered,  CODAR-type  HF  radar-derived  surface  currents  within  a  model  configuration  that  includes  tidal  forcing,  t  is  shown,  that 
assimilation  of  unfiltered  (with  tides)  surface-current  observations  into  the  model  with  tides  improves  the  sub-lidal  model  predictions  to  the  level  con  parable  with  the 
assimilation  of  filtered  data  into  the  non-tidal  model,  which  is  significant  with  respect  to  options  for  designing  real-time  nowcast  and  forecast  syst.ms.  finally,  the 
approach  is  extended  and  evaluated  for  the  direct  assimilation  of  HF  radar-derived  radial  velocity  components.  Direct  assimilation  of  radial  currents  t  Uends  the  range 
of  influence  of  the  data  into  regions  covered  by  only  one  HF  radar  site. 


15.  SUBJECT  TERMS 

Ocean  data  assimilation  and  reanalysis;  Coastal  processes;  Currents;  Ocean  observing  systems;  Upwelling 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

18.  NUMBER 
OF 

PAGES 

13 

19a.  NAME  OF  RESPONSIBLE  PERSON 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

Igor  Shulman 

Unclassified 

Unclassified 

Unclassified 

UL 

19b  TELEPHONE  NUMBER  (Include  area  co  e) 
228-688-5646 

Standard  Form  2  )8  (Rev.  8/98) 


Prescribed  by  ANSI  Si  Z39  18 


Deep-Sea  Research  11  56  (2009)  149-160 


Contents  lists  available  at  ScienceDirect 

Deep-Sea  Research  II 

journal  homepage:  www.elsevier.com/locate/dsr2 


KEF-SU  RESEMICT 
FACT! 


Assimilation  of  HF  radar-derived  radials  and  total  currents 
in  the  Monterey  Bay  area 

Igor  Shulman  a-*,  Jeffrey  D.  Paduan  b 

J  Naval  Research  Laboratory.  Code  7331.  Building  7009.  Room  A 144.  Stennis  Space  Center.  MS  39S29-S004,  USA 
b  Naval  Postgraduate  School.  Monterey.  USA 


ARTICLE  INFO 


ABSTRACT 


Article  history: 

Accepted  16  August  2008 
Available  online  19  September  200B 

Keywords: 

Ocean  data  assimilation  and  reanalysis 

Coastal  processes 

Currents 

Ocean  observing  systems 
Upwelling 


Impact  of  HF  radar  surface-current  assimilation  on  ocean  circulation  model  predictions  in  the  Monterey 
Bay  area  is  studied  and  evaluated  during  the  time  frame  of  the  Autonomous  Ocean  Sa  npling  Network 
(AOSN-il)  experiment  (August-September  2003).  In  the  first  instance,  a  previously  d(  scribed  method 
for  assimilation  of  surface  current  data  is  applied  to  33-h  low-pass-filtered  data  and  a  r  )n-tidal  version 
of  the  circulation  model.  It  is  demonstrated  that  assimilation  of  surface  velocity  d  ta  significantly 
improves  the  surface  and  subsurface  correlation  of  model  currents  with  moored  current  observations. 
These  results  from  the  AOSN-II  period  illustrate  that  surface-current  assimilation  is  beneficial  even  in 
cases  for  which  very  high-resolution  (3  km)  atmospheric  forcing  is  utilized.  The  assimik  i  ion  approach  is 
also  tested  with  hourly,  unfiltered,  CODAR-type  HF  radar-derived  surface  currents  vithin  a  model 
configuration  that  includes  tidal  forcing.  It  is  shown,  that  assimilation  of  unfiltered  (wil  i  tides)  surface- 
current  observations  into  the  model  with  tides  improves  the  sub-tidal  model  predict  ons  to  the  level 
comparable  with  the  assimilation  of  filtered  data  into  the  non-tidal  model,  which  is  ;ignificant  with 
respect  to  options  for  designing  real-time  nowcast  and  forecast  systems  Finally,  le  approach  is 
extended  and  evaluated  for  the  direct  assimilation  of  HF  radar-derived  radial  velocity  c  mnponents.  The 
model  runs  that  included  assimilation  of  radials  from  at  least  two  HF  radar  sil  *s  show  better 
correlations  with  observations  than  the  non-assimilative  run,  especially  those  runs  that  included  radials 
from  the  Santa  Cruz  site.  Directions  of  radials  for  that  site  coincide  with  the  directk  is  of  dominant 
southward  flow  during  upwelling  events  and  the  northward  flow  during  relaxation  events.  Direct 
assimilation  of  radial  currents  extends  the  range  of  influence  of  the  data  into  regions  overed  by  only 
one  HF  radar  site. 
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1.  Introduction 

During  the  last  decade  high-frequency  (HF)  radar  systems  have 
been  installed  operationally  throughout  the  world.  Assimilation  of 
HF  radar  surface  currents  into  oceanic  models  has  been  a  subject 
of  a  number  of  studies  (Lewis  et  al.,  1998,  Breivick  and  Saetra, 
2001 ;  Oke  et  al.,  2002;  Kurapov  et  al.,  2003;  Paduan  and  Shulman, 
2004;  Wilkin  et  al.,  2005)  Surface-current  data  assimilation 
experiments  based  on  high-frequency  radar  observations  in 
summer  1999  and  2000  were  described  in  Paduan  and  Shulman 
(2004).  In  that  study,  low-pass-filtered  surface  currents  were 
assimilated  into  a  non-tidal  circulation  model  of  Monterey  Bay 
based  on  a  nested  implementation  of  the  Princeton  Ocean  Model 
(POM).  That  model  was  forced  with  either  the  91-km-resolution 
winds  from  the  Navy’s  Global  Atmospheric  Prediction  System 
(NOGAPS;  Rosmond  et  al.,  2002)  or  9-km  resolution  winds  from 
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the  Navy’s  Coupled  Ocean/Atmosphere  Mesosc  le  Prediction 
System  (COAMPS™;  Hodur  et  al.,  2002).  The  evalu  ated  assimila 
tion  scheme  consists  of  two  steps:  the  physical-space  statistical 
analysis  system  (PSAS)  is  used  to  derive  correction  -  to  the  model 
surface  velocity  based  on  comparisons  with  ob;  ?rved  surface 
currents.  Then  corrections  are  projected  downwai  1  through  the 
frictional  boundary  layer  assuming  that  the  mode  -data  velocity 
differences  at  the  surface  represent  the  top  of  a  onstant  eddy 
viscosity  Ekman  boundary  layer  (see  Paduan  and  S  lulman,  2004, 
and  Section  3  below).  The  underlying  hypothesis  in  his  procedure 
is  that  inadequate  wind-stress  forcing  can  be  par  ally  compen 
sated  by  adjusting  model  currents  toward  the  ob  erved  surface 
currents.  It  was  shown  that  assimilation  of  CODAR  type  HF  radar 
data  improved  model  simulations  at  mooring  local  ions  down  to 
120m  (which  was  well  below  the  depths  directly  influenced  by 
the  Ekman-layer-assimilation  procedure;  Paduan  ind  Shulman, 
2004). 

The  present  study  represents  a  follow-on  to  the  work  of 
Paduan  and  Shulman  (2004)  that  takes  advantaj  ?  of  the  data 
collected  around  Monterey  Bay  as  part  of  the  Autonomous  Ocean 


150 


/.  Shulman.  J.D  Paduan  /  Deep-Sea  Research  ll  S6  (2009)  149  160 


Sampling  Network  Experiment  (AOSN-ll)  in  August-September 
2003.  Also  important  is  the  initiation  of  a  high-resolution  (31<m) 
COAMPS™  atmospheric  model  nest  covering  the  central  California 
region  that  was  first  put  in  place  during  AOSN-ll  (Doyle  et  aL, 
2008).  In  this  study,  we  address  the  following  issues: 

1  Impact  of  assimilation  of  low-pass-filtered  HF  radar  surface 
currents  on  model  predictions  during  AOSN-ll  time  frame. 

2.  Impact  of  assimilation  of  unfiltered  HF  radar  surface  currents 
on  model  predictions  during  AOSN-ll  time  frame.  The  data 
assimilation  approach  of  Paduan  and  Shulman  (2004)  was 
designed  for  correcting  wind-driven,  sub-tidal  currents.  For 
this  reason,  the  33-h  low-pass-filtered  CODAR  data  were 
assimilated  into  the  circulation  model.  Oke  et  al.  (2002),  for 
example,  used  a  rather  more  computationally  expensive  pre¬ 
filtering  in  their  data  assimilation  scheme  to  minimize  shocks 
in  the  model  responses.  We  investigate  whether  the  assimila 
tion  of  unfiltered  data  improves  or  degrades  sub-tidal  model 
predictions  in  comparison  to  the  assimilation  of  pre-filtered 
observations. 

3.  Impact  of  assimilation  of  HF  radar  radials  on  model  predictions 
during  AOSN-ll  time  frame.  Benefits  of  direct  radial  velocity 
assimilation  include  the  ability  to  avoid  the  total  vector 
combination  step  and  errors  associated  with  geometric  dilu 
tion  of  precision  effects.  The  area  of  data  influence  can  be 
extended  by  including  some  information  in  regions  covered  by 
just  one  HF  radar  site.  Direct  assimilation  of  radial  velocity 
components  also  expands  the  possibility  to  assimilate  HF  radar 
data  from  ships  and  petroleum  platforms.  Formulation  of  the 
assimilation  problem  in  terms  of  radial  velocity  components 
does,  however,  greatly  expand  the  apparent  number  of 
observations  that  must  be  dealt  with  during  each  assimilation 
time  step, 

2,  Observation  and  model  descriptions 

2.1.  HF  radar  network 

Surface-current  observations  used  in  this  study  were  derived 
from  a  network  of  five  SeaSonde-type  HF  radar  instruments 
deployed  in  the  region  of  Monterey  Bay  (Fig.  1).  Those  instru- 
ments,  commonly  referred  to  as  CODAR  type  HF  radar  systems, 
exploit  information  in  the  radiowave  backscatter  from  the  ocean 
surface  to  infer  movement  of  the  near  surface  water.  Electro¬ 
magnetic  waves  in  the  HF  band  (approximately  3*30  MHz)  exhibit 
Bragg-resonant  reflections  from  wind  driven  gravity  waves  on  the 
ocean  surface  whose  physical  wavelength  is  precisely  1/2  the 
wavelength  of  the  transmitted  radiowave.  During  the  AOSN-ll 
period  in  August-September  2003,  four  SeaSonde  systems  were 
operating  at  frequencies  near  12.5  MHz  and  one  system  (in  Moss 
Landing)  was  operating  at  25.4  MHz,  which  meant  that  the  Bragg- 
resonant  scatter  from  the  sea  surface  was  due  to  gravity  waves 
whose  wavelengths  were,  approximately.  12  and  6  m,  respectively. 

Several  studies  have  investigated  the  performance  of  the 
Monterey  Bay  HF  radar  network  by  comparing  the  radar-derived 
currents  with  in  situ  velocity  observations  and  by  comparing 
radar-to-radar  velocity  estimates  on  the  over-water  baselines 
between  radar  sites  (e.g.,  Paduan  and  Rosenfeld,  1996;  Paduan 
et  al.,  2006).  Consistent  uncertainty  values  emerge  in  the  range  of 
10-15  cm/s  for  the  remotely  estimated  velocities.  In  addition  to 
those  performance  measures.  Paduan  and  Shulman  (2004) 
described  monthly  tabulations  of  cross  shore  and  along  shore 
velocity  decorrelation  scales  based  on  earlier  computations  from 
the  Monterey  domain.  These  same  uncertainty  and  decorrelation 
values  are  used  in  this  study. 


123.4  122.6  12I.8W 


Fig.  1.  The  ICON  model  domain  with  local  bathymetry  and  the  loca  ions  of  coastal 
HF  radar  sites  (tnangles)  and  offshore  moorings  (Ml  and  M2). 


Also  relevant  to  this  study  are  the  basic  descript  ons  of  data 
availability  from  the  HF  radar  network.  Each  individt  il  SeaSonde 
instrument  provides  a  distribution  of  so-called  "rad  al”  velocity 
observations  each  hour  on  a  polar  coordinate  grid  cen  ered  on  the 
radar  site  Independent  estimates  of  the  speed  o  the  water 
approaching  or  receding  from  the  radar  site  are  provic  ed  at  scales 
of  3  km  in  range  (1.5  km  for  the  25.4  MHz  system  and  5  in 
azimuth.  Each  hour’s  spatial  set  of  radial  velocity  esti  nates  is  not 
necessarily  filled  in.  This  is  due  to  limitations  of  re  direction 
finding  algorithm  used  with  a  compact  HF  radar  sys  em  such  as 
the  SeaSonde  (see,  for  example,  Barrick  and  Lipa,  1997  Laws  et  al., 
2000;  de  Paolo  and  Terrill,  2007,  Toh,  2005).  The  cumulative 
radial  velocity  coverages  are  shown  in  Fig.  2  for  each  i  adar  site.  In 
the  figure,  the  value  at  each  grid  location  depicts  the  p  -rcentage  of 
the  total  possible  hourly  observations  obtained  at  t  iat  location 
during  the  analysis  period.  From  the  figure,  it  can  be  s  -en  that  the 
offshore  range  for  the  12.5-MHz  systems  was  betwe  n  about  50 
and  60km,  while  the  range  for  the  25.4-MHz  systen  was  about 
40km.  Vector  current  estimates  require  overlapping  radial 
observations  from  two  or  more  HF  radar  sites,  whin  results  in 
more  limited  coverage.  Vector  currents  were  estimated  on  a 
Cartesian  grid  with  a  horizontal  resolution  of  3  km  b\  computing 
the  best  fit  vector  velocity  components  using  all  radial  velocity 
observations  within  a  radius  of  3  km  for  each  grid  point  each  hour 
(hence,  neighboring  vector  current  results  are  not  completely 
independent).  The  percent  coverages  by  grid  locatio  i  for  (total) 
vector  currents  are  also  shown  in  Fig.  2. 

Finally,  it  is  important  to  point  out  what  are  the  <  pproximate 
depths  of  the  HF  radar-derived  current  estimates.  Bee  luse  the  HF 
radar  measurement  depends  on  the  ocean  currents  npeding  or 
assisting  the  Bragg-resonant  gravity  waves,  the  depth  >r  thickness 
of  the  relevant  ocean  currents  depends  on  the  penetration  depths 
of  the  resonant  wave's  particle  motions.  The  weight  -d  averaged 
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Fig.  2.  Radial  footprints  and  percent  temporal  coverage  for  HF  radar  sites:  SCRZ  (A;  Santa  Cruz),  MLNG  (B;  Moss  Landing),  PP1N  (C;  Pt.  Pinos).  NPGS  (D;  N  ival  Postgraduate 
School),  and  PSUR  (F;  Pt.  Sur).  The  footprint  and  temporal  coverage  for  the  vector  analysis  grid  are  shown  along  with  the  locations  of  moorings  Ml  an  l  M2  (E). 


depths  of  the  measured  currents  are  on  the  order  of  0.8  and  0.5  m 
for  HF  radar  measurements  at  12.5  and  25.4  MHz,  respectively 
(Stewart  and  Joy,  1974).  Note:  in  this  study  we  have  ignored  subtle 
differences  expected  for  currents  measured  at  these  very  near 
surface  depths.  In  other  cases,  investigators  have  attempted  to  use 
simultaneous  measurements  over  a  range  of  frequencies  to  look 
for  statistical  differences  in  the  measured  velocities  (e.g.,  Teague 
et  al .  2001 ;  Meadows,  2002). 


2.2.  Moored  observations 

For  this  study  (August  6  September  6  of  2003).  the  model 
current  predictions  with  and  without  assimilation  of  HF  radar 
data  were  compared  to  currents  measured  by  a  75-kHz  RD 
Instruments  Acoustic  Doppler  Current  Profiler  (ADCP)  mounted 
in  a  downward-looking  configuration  on  the  Monterey  Bay 


Aquarium  Research  Institute’s  (MBARI)  surface  moorings:  at 
122.02  W,  36.74  N,  designated  Ml  and  at  122.4  )  W,  36.67  N, 
designated  M2  (see  Fig.  1).  For  both  moorings  the  ADCPs  were  set 
up  to  measure  currents  in  60  bins  (8  m  bins  with  ie  first  bin  at 
16  m  depth).  Complex  correlation  and  mean  speed  e  rrors  between 
the  observed  and  model-predicted  currents  were  us  d  for  compar 
isons.  The  magnitude  p  and  the  angular  displacement  0  of  the 
complex  correlation  coefficient  between  the  ADCP  ,ind  the  model 
currents  for  a  particular  depth  were  estimated  usinj  the  approach 
outlined  in  Kundu  (1976).  The  magnitude  is  estimated  as: 

p  =  v  Re2  +  Im2.  (1) 

where 

\JYLMK12 + +  (vf1)2) 
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lm  E,(u?vr  +  y?uH 

v/Er«u?)2  +  (v?)2)Er«“!")2  +  <v!")2) 

The  corresponding  angular  displacement  0 ,  which  is  also  called 
the  phase  angle,  is  computed  according  to: 


0  =  tan'1 


Er(nfv™  -  vfu?) 
£r(u?u;"  +  v?v™) 


(2) 


where  utm,  vfm  are  the  demeaned  east-west  and  north-south 
model  velocity  components,  respectively,  and  ut°.  vr°  are  the 
demeaned  east-west  and  north-south  observed  velocity  compo¬ 
nents,  respectively.  The  angular  displacement  9  gives  the  average 
counterclockwise  angle  difference  between  model  and  observed 
velocity  vectors. 

In  addition  to  the  complex  correlation,  we  compute  a  measure 
of  error  based  on  the  average  of  the  vector  speed  difference 
according  to: 


RMSE  = 


(a^-un2+(v?-^)2),/2 

N 


(3) 


where  Utm,  and  Ut°yt°  are  the  total  velocity  components,  i.e., 
not  demeaned.  N  is  the  number  of  temporal  samples. 


2.3.  The  Monterey  Bay  model 

The  hydrodynamic  model  used  in  this  study  was  described  in 
detail  in  Shulman  et  al.  (2002),  Paduan  and  Shulman  (2004), 
Ramp  et  al.  (2005),  and  Rosenfeld  et  al.  (2008).  It  is  called  the 
ICON  model,  which  is  a  fine  resolution  ocean  model  covering 
Monterey  Bay  and  the  adjacent  coastal  areas  (see  model  domain 
on  Fig.  1 ).  The  horizontal  resolution  ranges  from  1-4  km  and  there 
are  30  vertical  sigma  levels. 

In  this  study,  the  ICON  model  was  forced  with  3-km  resolution 
wind  stresses  and  heat  fluxes  from  the  COAMPS  predictions.  Tidal 
forcing  was  introduced  into  the  model  through  the  specification  of 
the  open  boundary  conditions  using  tidal  constants  interpolated 
from  the  Oregon  State  University  Tidal  Solution  (Egbert  and 
Erofeeva,  2002)  for  the  US  West  Coast.  Eight  tidal  constituents 
(M2,  S2,  N2,  K2,  Kl,  01,  PI,  Q1 )  were  included. 

On  open  boundaries  the  ICON  model  is  coupled  to  the  Pacific 
West  Coast  (PWC)  model  (Haidvogel  et  al..  2000;  Rochford  and 
Shulman,  2000)  using  the  following  Flather  condition  (Rochford 
and  Shulman,  2000): 

Th  =  u°  +(g/H)'  (4) 

with 

n°  =  '/pwc  +  9  tide 

U°  =  u^T  +  u-^,  (5) 

where  rj  is  the  model  sea  surface  elevation  on  the  open  boundary, 
un  is  the  model  vertically  averaged  outward  normal  component  of 
velocity  on  the  open  boundary,  >/pwc  and  are  a  sea  surface 
elevation  and  vertically  averaged  velocity  from  larger  scale  PWC 
model  simulations  (no  tides).  Respectively,  r]tidc  and  are  a 
tidal  sea  surface  and  vertically  averaged  tidal  velocity  predicted  by 
using  eight  tidal  constants  from  the  OSU  model  (Egbert  and 
Erofeeva,  2002).  Details  about  the  tides  implementation,  along 
with  comparisons  of  vertical  and  horizontal  patterns  against  those 
derived  from  moorings  and  HF  radar  are  presented  in  Rosenfeld 
et  al.  (2008). 


3.  HF  radar  data  assimilation  approach 

According  to  Paduan  and  Shulman  (2004),  the  analysis 
(updated)  field  of  the  model  surface  currents  is  derived  from 


(bold  letters  are  vectors); 

Ua3  =  U' +K(U^-HUf,).  (6) 

where  Ua*  are  the  analyzed  surface  currents,  u/  are  the  model 
forecast  surface  currents,  U^°  are  HF  radar-derived  surface 
currents,  a  indicates  model  grid  locations,  and  p  indicates 
observations  locations.  H  is  the  interpolation  operator  of  the 
model  velocity  from  the  model  grid  a  to  the  observation  locations 
P  and  the  matrix  K  is  the  Kalman  gain,  which  depends  on  the 
forecast  error  covariance  matrix  Pf  and  the  observation  error 
covariance  matrix,  R: 

K  =  PfHT(HPrHT  +  Rr'.  (7) 

According  to  Eq.  (6),  corrections  to  the  surface  ve  3city  = 
Uj  —  Uj  =  (<5uj  ,<5va)  were  applied  to  the  surface  layer  of  the 
model.  These  surface  velocity  corrections  were  assumed  to  be 
related  to  errors  in  the  model  wind  forcing  and  were  projected 
downward  based  on  Ekman  theory  (for  more  details  see  Paduan 
and  Shulman,  2004).  F*  and  R  in  Eq.  (7)  are  describe  1  in  Paduan 
and  Shulman  (2004),  and  estimated  from  observed  HF  radar 
covariance  scales 

The  above  data  assimilation  approach  is  designed  for  correct¬ 
ing  wind-driven,  sub-tidal  currents.  All  experiment  in  Paduan 
and  Shulman  (2004)  were  conducted  using  33-h  low  tass-filtered 
CODAR  data  assimilated  into  the  circulation  model  v  ithout  tidal 
forcing.  With  the  implementation  of  tidal  forcii  g  into  the 
Monterey  Bay  model,  assimilation  of  unfiltered  surf  ce  currents 
became  more  consistent  with  the  model  forcing  fund  ions.  In  the 
case  of  the  model  with  tides,  errors  in  tidal  predicts  ns  will  also 
contribute  to  the  misfits  between  observed  and  modiM  predicted 
surface  currents  (Rosenfeld  et  al.,  2008).  In  this  study  the  same  P* 
and  R  were  used  for  assimilation  of  filtered  and  unfilt  :kred  surface 
currents  (the  objective  is  to  compare  sub  tidal  mode  predictions 
with  assimilation  of  filtered  and  unfiltered  data).  The  question 
here  is  whether  the  assimilation  of  unfiltered  data  mproves  or 
degrades  sub-tidal  model  predictions.  Being  fully  c  ware  about 
possible  drawbacks  of  the  approach,  answers  to  the  above 
question  are  addressed  in  Section  4.1. 

As  was  stated  in  the  Introduction,  another  obje  tive  of  the 
paper  is  development  and  testing  of  an  approach  f<  r  direct  HF 
radar  radial  current  assimilation.  In  that  case,  obse  vations  are 
projections  of  unknown  observed  vectors  UjJ  =  {U£,  V  }}  on  radial 
directions.  The  result  is: 

U^  =  (U“eR)  =  U“exR  +  V“eyR  (8) 

where  eR  is  a  unit  vector  along  the  radial,  and  U$,p  is  ,1  projection 
of  vector  U^°  on  the  radial  direction.  In  (6)  the  operaior  H  is  just 
the  interpolation  operator  of  the  model  velocity  to  the  locations  of 
observations.  In  the  case  of  radials  assimilation,  we  ii  troduce  the 
operator  HR,  which  does  the  interpolation  of  the  me  del  velocity 
from  the  model  grid  oc  to  the  observations  location  p  and  the 
projection  of  velocity  onto  the  direction  of  the  co.  responding 
radial.  Therefore,  the  operator  HR  both  interpolates  model  velocity 
and  transforms  vectors  to  scalars.  The  adjoint  operator  Hr 
transforms  the  scalar  to  a  vector  at  the  observatioral  locations 
(hence  the  use  of  bold  case  for  Hr)  and  interpolate  the  vector 
back  to  the  model  grid  locations: 

(HRu;.u^)  =  (ul.HTRu^. 

where 

hJur  =  U£eR 

In  this  case,  instead  of  (6)  we  have: 


ua  =  u'  +  KR(U7  -  Hr ufa>. 


(9) 


/.  Shufman.  J  D.  Paduan  /  Deep-Sea  Research  II  56  (2009)  149-160 


153 


where,  instead  of  (7),  we  have  the  following  expression  for 
Kalman  gain  KR. 

Kr  =  PfHj(HRPfHT  +  RPr\  (10) 

where  Rp  is  the  observational  error  covariance  for  radials.  In 
practice,  we  solve  first  a  linear  system  with  unknown  quantity  q 
such  that: 

(HRP‘Hl  +  Rp)q  =  Uy  -HrU[  (11) 


and  then  the  analyzed  state  U/  is  obtained  from  the  equation: 
u;  =  UfJ+PfHTRq.  (12) 


4.  Data  assimilation  experiments 

To  address  the  questions  that  we  posed  with  regard  to  pre 
filtering,  wind  resolution,  and  radial  current  assimilation,  we 
conducted  a  series  of  model  simulations.  All  analyses  were 
conducted  for  the  period  6  August-September  2003. 


122.3  122.2  122  1  122.0  121.9W  122.3  122.2  122.1  122.0  121.9  121.8 


Fig.  3.  CODAR  surface  currents  averaged  over:  (A)  entire  time  frame;  (B)  first  upwelling  (August  6-19);  (C)  first  brief  relaxation  (August  20-22);  (D  second  extended 
upwelling  (August  23-31)  and  (E)  second  brief  relaxation  (5eptember  1  3).  The  symbol  in  each  panel  denotes  the  location  of  the  Ml  mooring  (•). 
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As  presented  in  Ramp  et  al.  (2008),  there  were  two  extended 
upwelling  events  (August  6-19,  and  August  23-31 ),  two  very  brief 
relaxation  events  (August  20-22,  and  September  1-3),  and  a  weak 
upwelling  (September  4-6)  event.  Fig.  3  shows  CODAR  surface 
currents  averaged  over  the  entire  time  frame  and  each  of  these 
identified  events.  During  extended  upwelling  events.  Fig.  3B  and  D 
indicate  development  of  a  strong  southward-flowing  jet  along  the 


Table  1 

Model  and  data  configurations  for  various  runs  presented  in  this  study 


Run  # 

Surface  current  input  data 

Tidal  forcing 
Included  in 
modeP 

Assimilation  Assimilation  Sources'* 

of  low-pass-  of  unfiltered? 

filtered? 

1 


2 

Yes 

- 

Vector 

- 

3 

- 

- 

- 

Yes 

4 

- 

Yes 

Vector 

Yes 

5 

Yes 

- 

SCRZ,  MLNG. 
PPIN,  NPGS 

- 

6 

Yes 

- 

SCRZ,  MLNG 

- 

7 

Yes 

- 

MLNG,  PPIN 

- 

8 

Yes 

- 

SCRZ,  PPIN 

- 

9 

Yes 

- 

SCRZ.  PSUR 

- 

10 

Yes 

- 

MLNG 

- 

11 

Yes 

- 

SCRZ 

- 

12 

Yes 

- 

PPIN 

- 

13 

Yes 

PSUR 

- 

14 

Yes 

SCRZ,  MLNG. 
PPIN.  NPGS. 
PSUR 

J  Vector  data  was  formed  by  combining  radial  components  from  four  sites: 
SCR2.  Ml.NG.  PPIN.  NPGS. 


entrance  to  the  Bay.  Note,  that  CODAR  mean  currents  over  the 
entire  month  are  also  southward  near  the  Ml  mooring  (Fig.  3A). 
During  the  brief  relaxation  periods,  when  upwelling-favorable 
winds  weaken,  the  structures  are  more  complicated  yet  there  is 
clear  indication  of  the  development  of  along-shore  northward 
flow  along  the  entrance  to  the  Bay  (Fig.  3C  and  E). 

The  distinct  attributes  of  each  model  run  are  given  in  Table  1. 
All  simulations  were  conducted  using  the  same  initial  conditions 
and  were  forced  with  wind  stress  and  heat  fluxes  from  the  3-km~ 
resolution  COAMPS™  atmospheric  model  analyses. 

4  1.  Assimilation  of  HF  radar  surface  currents 

Run  1,  which  is  without  assimilation  of  surface-cu  rent  data  in 
the  non-tidal  model,  is  the  control  run  for  comparisor  s.  Run  2  is  a 
clone  of  Run  1  but  filtered  CODAR  surface  cu  rents  were 
assimilated.  Run  3  is  a  clone  of  Run  1  but  with  tic  es  included 
into  the  model  forcing,  and  Run  4  is  a  clone  of  Rur  3  but  with 
assimilation  of  unfiltered  CODAR  surface  currents  int<  the  model 

The  33-h  low-pass-filtered  observed  and  mod  4-predicted 
currents  at  48  m  depth  are  shown  in  Fig.  4.  A  qualit  tive  review 
shows  that  Runs  2  and  4  with  assimilation  of  HF  ndar  surface 
currents  reproduced  much  better  observed  subsurface  structure  of 
currents  at  the  Ml  mooring. 

Complex  correlations  and  angular  displacements  l  etween  the 
model  predicted  and  observed  currents  are  shown  at  mooring 
locations  Ml  and  M2  in  Figs.  5  and  6,  respectively. 

Comparisons  of  complex  correlations  and  angul  ir  displace¬ 
ments  for  Run  1  and  Run  2  show  that  assimilatio  of  filtered 
CODAR  surface  currents  into  the  non-tidal  circuh  tion  model 
significantly  improves  model  current  predictions.  /  :cording  to 
Fig.  5A,  on  average,  the  correlations  with  Ml  currenls  increased 
from  around  0.3-0.35  (Run  1)  to  around  0.7  (Run  2  Also,  with 


AUG  2003 

06  10  14  IK  22  26  30  03 


Fig.  4.  33-h  low-pass- filtered  observed  and  model-predicted  currents  at  48  m  depth  plotted  every  4h  for  the  Ml  mooring  location. 
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Fig.  5.  Complex  correlations  (A)  and  angular  displacements  (B)  between  model-simulated  and  observed  currents  at  mooring  Ml.  The  legend  refers  to  me  del  runs  described 
in  Table  1  and  the  notation  "filtered’’  indicates  that  both  model-simulated  and  observed  currents  were  low  pass  filtered  (33-h  half-power  point)  prior  t<  calculation  of  the 
complex  correlations. 


0  0  1  0.2  0.3  0.4  0.5  0.6  0.7  0  8 
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Fig.  6.  As  for  Fig.  5  but  at  mooring  M2. 
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assimilation  the  angle  between  model  currents  and  observations 
(Fig.  5B)  became  smaller  in  the  upper  40  m. 

As  we  stated  above,  there  were  five  upwelling/relaxation 
events  during  the  considered  time  frame.  With  five  degrees  of 
freedom,  a  correlation  value  of  0.67  is  significant  with  95% 
confidence  (see  for  example,  Table  XI  of  Hogg  and  Tanis  (1993)). 
With  this  interpretation,  correlations  are  significant  down  to  70  m 
depth  for  Runs  2  and  4  (filtered).  However,  we  can  suppose  that 
the  number  of  degrees  of  freedom  is  larger  than  five  based  on  the 
shorter  term  fluctuations  in  the  data  set.  There  is  strong  diurnal 
variability  in  the  atmospheric  conditions  in  the  area;  correlations 
were  estimated  over  a  31 -day  period  using  hourly  data  pass 
through  a  33-h  low-pass  filter.  This  gives  around  31*24/33~22 
samples.  Therefore,  the  actual  number  of  degrees  of  freedom  is 
between  5  and  22.  With  22  degrees  of  freedom  a  correlation  of 
0.34  is  significant  at  the  95%  confidence  level.  With  that 
interpretation,  correlations  for  Runs  2  and  4  (filtered)  are 
significant  for  all  depths  shown  on  Fig.  5,  and  down  to  100  m  for 
Run  4  (unfiltered).  Also,  correlations  for  Run  1  are  close  to  this 
significance  level  for  many  depths. 

Note,  that  the  value  of  angular  displacement  0  is  only 
meaningful  if  complex  correlation  p  is  significant. 

While  correlations  at  Ml  illustrate  the  impact  of  CODAR  data 
assimilation  on  the  model  predictions  inside  the  HF  radar 
footprint  (Fig.  1),  the  correlations  at  M2  illustrate  the  influence 
of  assimilation  on  model  predictions  outside  the  footprint. 
Comparison  of  Runs  1  and  2  in  Fig.  6  indicates  an  increase  in 
correlation  at  M2  for  the  run  with  assimilation.  This  corresponds 
with  results  presented  in  Paduan  and  Shulman  (2004),  where  it 
was  demonstrated  that  assimilation  of  HF  radar-derived  surface 
currents  improves  model  predictions  not  only  in  the  area  of  the  HF 
radar  coverage  but  also  outside  of  the  area  (for  distant  effect  of 
data  assimilation  see  also  Kurapov  et  al„  2005). 

On  Figs.  5  and  6,  curves  labeled  "Run  3  filtered”  and  '  Run  4 
filtered”  show  the  impact  of  assimilation  of  unfiltered  currents  on 
sub-tidal  predictions  of  the  model  with  tides.  There  is  a  significant 
improvement  in  predictions  of  sub-tidal  currents  with  the  model 
assimilating  unfiltered  data.  The  model  sub-tidal  currents  have 
similar  correlation  with  observations  in  the  case  of  assimilation  of 
filtered  data  into  the  non-tidal  model  (Run  2)  or  with  assimilation 
of  unfiltered  data  into  the  model  with  tides  (Run  4,  filtered).  Note 
that  at  mooring  M2,  correlation  for  "Run  4  filtered”  is  even  better 
than  for  Run  2. 

Correlation  results  are  well  supported  by  the  RMSE  values 
presented  on  Fig.  7  for  the  Ml  and  M2  locations.  There  is  a 
decrease  in  RMSE  in  sub  tidal  predictions  due  to  assimilation  of 
filtered  as  well  as  unfiltered  surface-current  data  into  the  model 
without  or  with  tides  in  comparison  to  the  non-assimilative  Run  1. 
This  result  is  encouraging  for  real  time  forecasting  systems  for 
which  pre-filtering  of  observations  is  expensive  and  sometimes 
impossible. 

The  inclusion  of  tidal  forcing  in  the  circulation  model  appears 
to  have  introduced  an  additional  source  of  error  with  respect  to 
the  mooring  observations.  This  may  not  be  surprising  given  that 
tidal  currents  in  the  region  vary  over  small  spatial  scales 
(Rosenfeld  et  al..  2008).  Location  or  phase  errors  in  the  modeled 
tidal  currents  can  be  expected  to  increase  model-observation 
velocity  differences.  Results  for  unfiltered  Runs  3  and  4  on  Figs. 
5-7  provides  an  estimate  of  the  impact  of  surface  current 
assimilation  on  total  model  current  predictions.  For  example, 
there  is  a  reduction  in  RMSE  and  a  much  better  correlation  for  Run 
4  with  the  surface  current  assimilation.  According  to  Figs.  5  and  6, 
assimilation  of  surface-current  data  in  Run  4  increased  correlation 
at  Ml  from  0.2  (Run  3)  to  0.7  at  the  surface  and  from  around  0.3  to 
0.5  on  average  in  the  subsurface;  at  M2  the  correlation  was 
increased  to  an  average  level  of  0,45  in  the  upper  100  m, 


Fig.  7.  RMSE,  computed  according  to  Eq.  (3),  between  the  model  simulated  and 
observed  currents  at  moorings  Ml  (left-hand  group)  and  M2  (rigl  hand  group). 
Dashed  lines  simply  highlight  the  two  sets  of  curves 

Surprisingly,  the  reductions  in  RMSE  due  to  assimilation  (Fig.  7) 
are  larger  at  M2,  which  is  just  outside  the  sur  ace-current 
coverage  area,  than  they  are  at  Ml,  which  is  in  the  c  rnter  of  the 
surface-current  coverage  (Fig.  2).  We  note  also  that  R  /1SE  for  Run 
3  with  tides  is  significantly  larger  than  that  for  the  un  without 
tides  (Run  1 ). 

To  insure  that  the  upper  ocean  results  reflect  the  full  impact  of 
surface  current  assimilation,  velocity  difference  stanstics  were 
investigated  for  a  subset  of  the  model  runs  using  11  available 
observation  depths.  Correlation  magnitude  and  phas  results  for 
the  control  run  (Run  1)  versus  the  assimilation  run  Run  2)  are 
shown  in  Figs.  8  and  9  for  the  Ml  and  M2  locations,  espectively. 
RMSE  values  for  both  mooring  locations  are  showi  in  Fig.  10. 
Based  on  the  correlation  results,  the  positive  effect  of  surface- 
current  assimilation  extends  to  about  140m  at  Ml  and  to  about 
200m  at  M2  Looking  at  the  RMSE  values  in  Fig.  10,  it  s  clear  that 
the  absolute  comparison  between  observed  anc  simulated 
currents  is  better  at  Ml  than  at  M2  for  all  depths.  \t  mooring 
Ml,  the  data  assimilative  run  has  smaller  values  of  RI\  SE  down  to 
about  270m  depth,  and  slightly  larger  values  belo  v  270m  in 
comparison  to  the  non-assimilative  run.  At  mooring  1 12,  the  data 
assimilative  run  has  smaller  values  of  RMSE  for  a  depths  in 
comparison  to  the  non-assimilative  run.  Overall,  it  appears  that 
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Fig.  8.  Complex  correlations  (solid)  and  angular  displacements  (dashed)  between 
model-simulated  and  observed  currents  at  mooring  M 1  for  Run  1  (bold )  and  Run  2 
(thin).  Both  model-simulated  and  observed  currents  were  low  pass  filtered  (33-h 
half-power  point)  prior  to  calculation  of  the  complex  correlations. 
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Fig.  9.  As  for  Fig.  8  but  at  mooring  M2. 

the  influence  of  surface-current  assimilation  is  limited  to  the 
upper  ocean  within  and  above  the  main  thermocline  consistent 
with  the  correction-driven  surface  divergence  mechanism  pro¬ 
posed  by  Paduan  and  Shulman  (2004),  while  the  depth  and 
magnitude  of  those  influences  can  vary  substantially  for  a  given 
domain. 


4.2.  Assimilation  of  HF  radar  radials 

Finally,  we  discuss  results  of  the  assimilation  of  radials  from 
the  five  HF  radar  sites  shown  on  Fig.  2  We  note  first  that  Run  2 
and  Run  4,  which  were  discussed  above,  assimilated  preprocessed 
vector  surface  currents  based  on  input  from  four  of  the  five  sites: 
SCRZ,  PPIN,  MLNG,  and  NPGS.  The  resulting  vector  grid  for  those 


Fig.  10.  RMSE,  computed  according  to  Eq.  (3),  between  the  mo  fel-simulated  and 
observed  currents  at  moonngs  Ml  (left-hand  group)  and  M2  ght-hand  group) 
for  Runl  (bold)  and  Run2  (thin).  Both  model-simulated  and  ubserved  currents 
were  low  pass  filtered  (33-h  half-power  point)  prior  to  calculat  )n  of  the  RMSE. 


data  is  shown  in  Fig.  2E.  Table  1  describes  the  alinbutes  of  all 
model  runs,  including  Run  5  through  Run  14  n  which  we 
assimilated  radial  surface  current  data  from  diffe  ent  combina¬ 
tions  of  the  available  HF  radar  sites.  In  these  radia  current  runs, 
filtered  observations  of  radials  currents  were  assimilated  into  the 
model  without  tides. 

Brief  descriptions  of  physical  conditions  in  the  a  ea  were  given 
in  Section  4.1  and  Figs.  3  and  4.  For  successful  i  ^simulation  of 
radial  current  data,  it  is  important  to  resolve  the  do  ninant  along¬ 
shore  component  of  the  observed  flow  during  upv  elling/relaxa 
tion  events,  especially  the  southward  flowing  jet  di  ring  extended 
upwelling  events.  Based  on  this,  and  the  footprint  ol  HF  radar  sites 
(Fig.  2),  it  is  clear  that  sites  MLNG,  NPGS,  an  1  PSUR  have 
limitations  in  resolving  the  flow  around  mooring  Ml.  The  MLNG 
radials  are  mostly  oriented  in  the  cross-shore  dire i  tion  and  are, 
therefore,  orthogonal  to  the  direction  of  the  domina  t  along-shore 
flow.  Mooring  Ml  is  located  close  to  the  boundari  *s  of  the  data 
footprints  for  the  NPGS  and  PSUR  sites.  Therefore,  :mly  the  SCRZ 
and  PPIN  footprints  cover  the  area  around  Ml  w  41.  While  the 
direction  of  the  SCRZ  radials  coincides  nicely  with  t  e  direction  of 
dominant  along-shore  flow  around  the  Ml  moorir  ;  (Fig.  3),  the 
PPIN  site  radials  veer  toward  a  cross-shore  directn  n  around  Ml 
and,  therefore,  are  expected  to  lose  information  about  the  major 
along-shore  flow  component.  Also,  the  SCRZ  foot  rint  provides 
coverage  of  currents  in  extended  areas  upstrea  ^  as  well  as 
downstream  of  Ml,  while  the  PPIN  site  provides  lin  ited  coverage 
of  the  area  to  the  south  Pt.  Pinos  where  northward  low  develops 
during  relaxation  events. 

Fig.  11  shows  complex  correlation  magnitude ;  for  various 
model  runs  (Table  1)  based  on  comparisons  v  ith  observed 
currents  at  Ml  For  reference,  results  for  Rui  1  (without 
assimilation)  and  Run  2  (with  assimilation  of  preprocessed 
surface  currents)  are  also  included.  Overall,  complex  correlation 
magnitudes  for  the  runs  with  assimilation  of  adial  surface 
currents  are  lower  than  those  for  Run  2  (with  a  simulation  of 
preprocessed  surface  currents)  at  the  Ml  location.  A  it  is  stated  in 
Section  2.1,  HF  radar  vectors  (before  being  assimilate,  in  the  Run  2) 
were  estimated  by  computing  the  best-fit  vt  ctor  velocity 
components  using  all  radial  observations.  In  this  case,  around 
mooring  Ml,  preprocessed  HF  radar  surface  curn  nts  represent 
better  observed  features  than  separate  radials  from  multiple  sites 
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Fig.  11.  Complex  correlations  between  model-simulated  and  observed  currents  at 
mooring  Ml.  The  legend  refers  to  model  runs  described  in  Table  1. 

(which  are  treated  as  independent  observations  in  the  data 
assimilation  scheme  of  Section  3).  At  the  same  time,  complex 
correlation  magnitudes  for  the  runs  with  assimilation  of  radials 
are  higher  than  those  for  non-assimilative  run  (Run  1)  at  the  Ml 
location.  Only  Run  10  (assimilation  of  only  MLNG  data)  and  Run 
12  (assimilation  of  only  PP1N  data)  have  complex  correlation 
magnitudes  that  are  essentially  the  same  as  those  for  the  non 
assimilative  run.  As  was  hypothesized  above,  radials  for  the  MLNG 
and  PP1N  sites  do  not  resolve  the  dominant  alongshore  flow  very 
well.  Overall,  model  runs  that  included  assimilation  of  radials 
from  at  least  two  HF  radar  sites  show  a  better  correlation  than  Run 
1,  especially  those  runs  that  included  radial  surface  current  data 
from  the  SCRZ  site.  This  is  in  accord  with  the  above  hypothesis 
that  the  directions  of  radials  for  the  SCRZ  site  coincide  with  the 
direction  of  dominant  southward  flow  during  upwelling  events 
and  the  northward  flow  during  relaxation  events. 

According  to  statistical  significance  discussions  in  Section  4  1, 
under  an  assumption  of  five  degrees  of  freedom  (five  upwelling/ 
relaxation  events)  correlations  are  significant  down  to  70  m  depth 
for  all  runs  with  at  least  two  sites  assimilated  (except  Run  7  with 
MLNG  and  PP1N  sites  assimilated).  Assuming  22  degrees  of 
freedom  (see  explanations  in  Section  4.1)  correlations  for  all  runs 
(Fig.  11 )  are  significant  at  least  to  100m  depth,  except  for  the  non- 
assimilative  run,  and  runs  assimilating  radials  from  the  MLNG 
(Run  10)  and  PP1N  (Run  12)  sites. 

Fig.  12  shows  complex  correlation  magnitudes  for  various 
model  runs  based  on  comparisons  with  observed  currents  at  M2. 
These  results  differ  from  those  centered  on  the  Ml  mooring  site. 


Assimilation  of  radials  from  only  one  site  SCRZ  (Ri  i  11)  show 
better  correlations  with  observed  currents  at  M2  th  in  assimila¬ 
tion  of  preprocessed  surface  currents  in  Run  2  (Fig.  12  One  of  the 
reasons  for  this  is  that  mooring  M2  is  located  outsit  e  of  the  HF 
radar  observational  footprint  for  vector  currents  (Fit  2E),  while 
the  footprint  of  the  SCRZ  site  radials  covers  the  M2  1  cation  (Fig. 
2A).  All  Runs  (Fig.  12)  that  include  radial  data  from  tv  o  sites,  one 
of  which  is  SCRZ,  also  perform  better  at  the  M2  locati  n  than  Run 
2  (for  example,  Runs  6  and  8). 

Note,  that  the  radial  footprint  of  the  PP1N  site  also  reaches  the 
M2  location.  However,  as  shown  above,  the  direction  t  <f  radials  for 
PP1N  (Fig.  2C)  do  not  resolve  the  dominant  along  shore  flow, 
which  results  in  low  correlations  with  observed  cur  ents  at  the 
Ml  mooring.  From  Fig.  12,  results  from  the  assimilatu  n  of  radials 
from  only  the  PP1N  site  (Run  12)  are  worse  in  comp^  ison  to  the 
assimilation  of  SCRZ  radials  (Run  11),  but  are  compa  able  to  the 
results  of  Run  2. 

The  PSUR  site  radials  also  reach  the  M2  mooring.  However, 
PSUR  (Fig.  2F)  has  a  coarser  resolution  footprint  aro  ind  the  M2 
mooring  and  the  site  is  located  far  away  in  compari  on  to  SCRZ 
and  PP1N  sites.  As  a  result,  correlations  for  Ru  l  13,  with 
assimilation  of  radials  from  only  the  PSUR  site,  are  lower  than 
those  for  runs  with  assimilation  of  radials  from  SO  l  and  PP1N 
(Runs  11  and  12).  However,  the  addition  of  radial  dat,  from  PSUR 
in  Run  14  (when  data  from  all  sites  were  assimilate!  )  improved 
model  predictions  down  to  70  m  depth  in  comparis<  n  to  Run  5 
(using  four  other  sites  except  PSUR).  Overall,  the  assimilation  of 
radial  surface  current  data  from  four  or  five  sites  in  Run  5  and  Run 
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14,  respectively,  produced  much  better  correlation  with  currents 
observed  at  M2  than  did  Run  2,  which  used  vector  current  data 
whose  observational  footprint  does  not  reach  the  M2  location. 

The  assimilation  of  radial  surface-current  data  from  just  the 
MLNG  site  (Run  10)  led  to  a  dramatic  decrease  in  the  complex 
correlation  magnitudes,  especially  for  mid  depths  at  M2,  com¬ 
pared  to  Run  1  with  no  data  assimilation.  This  is  due  to  the  fact 
that  the  MLNG  radial  footprint  does  not  reach  mooring  M2,  and,  as 
it  was  stated  above,  radial  directions  do  not  resolve  the  dominant 
along-shore  flow. 


5.  Conclusions 

Paduan  and  Shulman  (2004)  conducted  data  assimilation 
experiments  with  input  of  33-h  low-pass-filtered  surface  current 
observations  into  the  ICON  circulation  model  without  tides.  The 
model  wind  forcing  had  either  very  coarse  resolution  (around 
91  km)  from  the  Navy’s  NOGAPS  atmospheric  circulation  model  or 
relatively  fine  resolution  (around  9  km)  from  the  Navy’s 
COAMPS1  M  model.  In  those  cases,  the  assimilation  improved 
model  current  predictions  at  the  Ml  and  M2  mooring  locations  in 
Monterey  Bay.  and  it  was  hypothesized  that,  at  least,  part  of  those 
results  could  be  attributed  to  the  HF  radar-derived  surface-current 
fields  acting  to  correct  for  missing  structure  in  the  wind  fields.  In 
the  present  paper,  the  approach  developed  in  Paduan  and 
Shulman  (2004)  was  applied  to  the  assimilation  of  HF  radar- 
derived  surface  current  data  during  the  AOSN-II  experiment  time 
frame  in  August  and  September  of  2003  when  high-resolution 
(3  km)  COAMPS™  wind  forcing  was  available.  It  was  shown  that 
assimilation  of  surface-current  data  has  positive  impacts  on 
surface  and  subsurface  model  predictions  even  in  cases  for  which 
high-resolution  atmospheric  forcing  is  used. 

Another  aspect  of  the  present  study  involved  an  assessment  of 
the  role  played  by  high-frequency  velocity  fluctuations.  Currents 
in  the  Monterey  Bay  study  area  include  significant  contributions 
from  internal  tide  fluctuations  at  the  semi-diurnal  period  and  sea 
breeze  related  fluctuations  at  the  diurnal  period  (Rosenfeld  et  al.. 
2008).  It  has  been  an  open  question  as  to  whether  the  inclusion  of 
these  velocity  constituents  within  the  data  assimilation  procedure 
would  degrade  the  results  achieved  for  sub-tidal  period  currents. 
The  results  presented  here  show  that,  under  the  present 
assimilation  scheme,  sub-tidal  period  velocity  simulations  are 
similarly  improved  through  the  assimilation  of  either  low- 
pass-filtered  surface  currents  or  instantaneous  (hourly)  surface 
currents.  This  result  is  consistent  with  the  mechanism  proposed 
by  Paduan  and  Shulman  (2004)  to  explain  the  subsurface  impacts 
of  surface  current  assimilation  They  suggested  that  model  depths 
below  the  surface  assimilation  layer  are  influenced  when  the 
velocity  correction  field  at  the  surface  is  horizontally  divergent. 
The  additional  simulations  described  here  also  suggest  that  the 
subsurface  impacts  are  driven  by  those  divergences  with  time 
scales  that  are  longer  than  the  local  inertial  period,  i.e.  longer  than 
the  geostrophic  adjustment  time  scale. 

The  option  for  direct  assimilation  of  radial  surface  current 
observations  from  multiple  HF  radar  sites  also  was  demonstrated 
in  the  simulations  described  here.  Overall,  complex  correlation 
magnitudes  for  the  runs  with  assimilation  of  radial  surface 
currents  were  lower  than  those  for  the  run  with  assimilation  of 
preprocessed  surface  currents  at  mooring  Ml,  and  higher  at 
mooring  M2.  The  model  runs  that  included  assimilation  of  radials 
from  at  least  two  HF  radar  sites  show  a  better  correlation  than  the 
non-assimilative  run.  especially  those  runs  that  included  radials 
from  the  Santa  Cruz  site.  Directions  of  radials  for  that  site  coincide 
with  the  direction  of  dominant  southward  flow  during  upwelling 
events  and  the  northward  flow  during  relaxation  events.  This 


result  has  practical  implications  in  that  single-component 
information  from  the  radial  current  observations  can  be  included 
from  those  portions  of  the  model  domain  observed  by  just  one  HF 
radar  site.  In  real-time  applications,  this  approach  will  also  allow 
processing  to  continue  smoothly  even  when  the  available  surface 
current  coverage  changes  dramatically  due  to,  for  example, 
communication  failures.  At  the  same  time,  the  impact  of  single 
component  radial  data  from  individual  HF  radar  sites  was 
extremely  variable.  These  facts  point  to  the  need  to  investigate 
further  the  impact  of  radial  surface  current  assimilation,  including 
the  development  of  appropriate  weighting  schemes  based  on  the 
independence  and  quality  of  individual  radial  current  observa 
tions.  Such  figures  of  merit  are  being  developed  foi  other  types  of 
HF  radar  applications  (e.g.,  Lipa,  2003;  Kaplan  and  Lekien,  2007), 
which  should  be  useful  in  future  assimilation  procedures. 

In  summary,  we  find  this  recently  completed  set  of  case  studies 
to  represent  an  encouraging  step  toward  the  p  actical  use  of 
remotely  sensed  surface-current  mapping  data  tn  lata  assimilat¬ 
ing  coastal  ocean  models.  Although  the  impact  of ;  jrface-current 
data  is  limited  to  the  upper  ocean  above,  about,  1,  Om,  improve¬ 
ment  of  model  circulation  fields  in  that  depth  ranj.e  is  important 
because  it  includes  the  euphotic  zone  and  man  of  the  most 
critical  transport  processes.  The  results  are  eniouraging  also 
because  the  coastal  ocean  area  being  mapped  by  HF  radar 
instruments  is  growing  rapidly  making  those  instr  iments  one  of 
the  essential  components  of  the  nation’s  inte  grated  ocean 
observing  system. 
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